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22.1  Introduction 

The  combination  of  numerical  models  and  observational 
data  can  provide  a  unique  tool  for  studying  the  complex 
interactions  of  the  atmosphere,  the  ice,  and  the  ocean.  The 
formulation  of  numerical  ice  and  coupled  ice-ocean-atmO' 
sphere  models  is  based  on  our  knowledge  of  dynamic  and 
thermodynamic  principles  and  how  they  relate  to  ob=erved 
ice  conditions.  Field  experiments  such  as  the  Arctic  Ice 
Dynamics  Experiment  (AIDJEX)  [Pritchard,  1 980]  and  the 
Marginal  Ice  Zone  Experiment  (MIZEX)  [Journal  of  Geo- 
physical  Research  Oceans,  88fC5J,  92(C7),  96(C3)]  have 
provided  observational  data  from  which  the  basis  of  many 
of  the  formulations  for  ice  drift,  internal  ice  stresses,  heat 
exchange,  etc. ,  have  come.  Numerical  models,  on  the  other 
hand,  may  be  used  to  provide  information  on  ice  drift,  ice 
thickness,  and  ice  concentration  in  regions  where  observa¬ 
tions  are  scarce  or  missing.  In  addition,  numerical  models 
may  be  used  to  forecast  ice  conditions. 

Satellites  have  been  able  to  provide  observational  data 
over  larger  areas  and  for  longer  periods  of  time  than  from 
conventional  observations.  Since  1972,  passive  microwave 
data  from  polar  orbiting  satellites  have  provided  large-scale 
coverage  of Arctic  and  Antarctic  sea  ice  on  a  nearly  continu¬ 
ous  basis  at  resolutions  as  fine  as  25  km.  Visible  and 
infrared  imagery  provide  higher  resolution  data  than  the 
passive  microwave,  but  do  not  have  the  all-weather  capabil¬ 
ity  that  the  passive  microwave  data  have.  A  very  promising 
source  of  high-resolution,  all-weather  data  is  the  satellite- 
borne  Synthetic  Aperture  Radar  (SA.R).  Data  from  this 
instrument  have  recently  become  available  via  the  launch 
of  the  ERS-1  satellite  on  July  17, 1991. 

With  the  recent  availability  of  larger  amounts  of  satellite 
data  in  ice  covered  regions  and  the  continuous  development 
ot  more  advanced  numerical  models  for  these  regions,  the 
logical  next  step  is  merging  the  data  with  the  models.  This 
merger  will  most  likely  occur  through  expanded  use  of  data 
assimilation  into  the  models.  This  combination  should 
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provide  the  best  available  analysis  and  forecast  of  ice 
conditions,  as  well  as  improve  our  understanding  of  dy¬ 
namic  and  thermodynamic  interactions  in  the  ice.  In 
addition,  models  are  likely  to  be  used  more  extensively  to 
improve  data  processing  algorithms  and  to  help  interpret 
complex  responses  recorded  in  the  data. 

22.2  Numerical  Models 

22.2.1  Ice  Models — Dynamic,  Thermodynamic,  and 
Dynamic-Thermodynamic 

Numerical  models  of  sea  ice,  developed  over  approxi¬ 
mately  the  last  twenty  years,  may  be  broken  down  into 
three  categories;  dynamic,  thermodynamic,  and  dynamic- 
thermodynamic  ice  models. 

Dynamicice  models  use  the  various  stresses  on  the  ice  to 
define  the  motion  of  the  ice.  The  momentum  balance  used 
in  dynamic  ice  models  is 

m  ~^u  =  -m  f  k  XU  +  Ta  +  Tu,.- mg  gradH  +  F  (1) 
5  t 

where  m  is  the  ice  mass  per  unit  area  and  ims  the  ice  drift 
velocity.  The  terms  on  the  right  side  of  the  equation  are  the 
Coriolis  term  (where /"is  the  Coriolis  parameter),  the  wand 
stress  on  ice  the  ocean  stress  on  ice  r^,  the  acceleration 
due  to  sea  surface  tiltmg  grad//(where  H  is  the  sea_^urface 
dynamic  height),  and  the  internal  ice  stress  term  F. 

The  simplest  of  the  dynamic  models,  the  free  drift  model, 
uses  a  balance  between  top  and  bottom  surface  stresses  and 
the  Coriolis  force  todetermine  ice  motion.  Free  drift  is  often 
a  good  approximation  away  from  boundaries  and  under 
divergent  conditions.  Removal  of  the  last  term  in  Equation 
(1),  the  effect  of  internal  ice  stress,  results  in  the  free  drift 
balance.  Ice  drift  may  be  significantly  adjusted  in  both 
magnitude  and  direction  by  the  internal  ice  stress  that 
generally  acts  in  a  direction  opposite  to  the  resultant  of  the 
wind  stress  and  Coriolis  force.  Models  that  include  internal 
ice  stress  contain  a  constitutive  law  that  treats  ice  as  a 
viscous,  elastic,  viscous-plastic,  or  ela.stic-pla,stic  medium. 
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In  a  viscous  rheologj',  stress  can  only  be  sustained  through 
a  nonrecoverable  dissipation  of  energy  by  deformation.  A 
plastic  rheology'  allows  stress  to  be  sustained  through  a  lack 
of  deformation  or  elastic  deformation  in  which  the  energy  is 
recoverable.  A  linear  viscous  rheology  in  which  stress 
depends  on  linear  strain  rates  and  a  rigid  plastic  rheology 
in  which  the  stress  stale  is  either  dependent  on  the  magni¬ 
tude  of  the  strain  rates  or  indeterminate  are  two  special 
case  rheologies  often  used  to  describe  sea  ice.  For  a  more 
detailed  review  of  these  constitutive  laws,  see  Hiblcr  [1980b, 
1986]. 

An  intermediate  approach  with  considerable  promise  for 
climate  simulations  is  the  cavitating  fluid  approximation, 
which  differs  from  free  drift  by  allowing  nonzero  ice  pres¬ 
sure  under  converging  conditions,  but,  like  free  drift,  offers 
no  resistance  to  divergence  or  shear  [Fiato  and  Hibler, 
1989],  Inthecorresponding  numerical  procedure,  free-drift 
velocities  are  corrected  itera  u  vely  in  a  momentum-conserv¬ 
ing  way.  ; 

Thermodynamic  ice  models  take  into  account  the  interac¬ 
tions  of  both  the  atmospheric  and  oceanic  heat  fluxes  with 
the  ice  to  determine  ice  growth  and  decay.  Many  of  the 
formulations  of  sea  ice  thermodynamics  used  in  numerical 
ice  models  are  based  on  the  worl ;  ofMaykut  and  Untersteiner 
[1 971  ] .  Figure  22-1  depicts  the  one-dimensional  snow-ice— 
water  system  they  use.  The  heat  fluxes  from  the  atmo¬ 
sphere  and  ocean  are  passed  into  the  snow  and  ice  by 
conduction. 

At  the  snow-air  surface,  the  balance  of  heat  fluxes  is 
given  by 

(1  -a}Fr-!o+FL-scT''  +Fs+F^  +  Ko 


0ifr<  273  K 
(S  +  /i)ifr>273  K 


(2) 


where  a  is  the  albedo,  is  the  solar  (shortwave)  radiation, 

and  Iq  is  the  flux  of  shortwave  radiation  through  the  surface 
into  the  ice,  assumed  by  Maykut  and  Untersteiner  to  be  1 7% 
of  the  net  shortwave  radiation  at  the  surface.  The  term  F/, 
is  the  incoming  longwave  radiation  from  the  atmosphere 
(clouds)  and  eoT^  is  the  outgoing  longwave  radiation  based 
on  the  surface  temperature,  where  e  is  the  longwave  emis- 
sivity,  cr  is  the  Stefan-Boltzmann  constant,  and  T  is  the 
surface  temperature.  Then,  F,  is  the  sensible  heat  flux, 
is  the  latentheat  flux,  Tifis  the  thermal  conductivity,^  is  the 
latent  heat  of  fusion,  z  is  the  depth  within  the  ice-snow 
column,  and  d{S  +  h}/dt  is  the  ablation  rate  of  snow  and  ice 
where  S  is  the  thickness  of  the  snow  and  h  is  the  thickness 
of  the  ice.  The  subscript  0  generically  refers  to  the  upper 
surface.  At  the  snow-ice  interface,  the  balance  of  fluxes  is 


h  SNOW 


Fig.  22-1  The  one-dimensional  snow-ice-water  s>-.sU;m  used  by 
Maykut  and  Untersteiner  [19711, 


ar.,  ^  ar, 

As-—  =Ki- — 
dz 
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The  subscripts  s  and  i  stand  for  snow  and  ice,  respec¬ 
tively.  At  the  ico-occan  interface  the  balance  is 


=Fu-q  —  {S^h)  14) 

()z  dt 

where  F„.  is  the  oceanic  heat  flux. 

Thermal  conductivity,  specific  heat,  and  ice  densjty  are 
all  functions  ofsalinity  and  temperature.  Thi.s  dependency 
is  due  in  large  part  to  brine  pockets  in  the  ice.  Brine  pocket, s 
can  act  as  thermal  reservoirs  and  .slow  down  the  heating  or 
cooling  of  ice.  Since  brine  has  a  small  conductivity  and 
greaterspecific  heat  than  ice,  the  paramoter.s  also  vary  with 
temperature. 
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To  calculate  heat  conduction,  the  diffusion  equation  for 
temperature  is  given  by 


The  model  also  incorporates  the  equations  of  continuity 
for  ice  concentration  and  thickness  defined  as 


(pc)i 


El 


-Ki 


■k^  lo  exp 


(5) 


El 

dt 


d(uh\ 

dx 


d(vh]  . 

— —  +  A/,  +  I  Difiusionj 

dy 


(7 1 


where  k  is  the  bulk  extinction  coefficient  and  c  is  the  heat 
capacity.  In  snow,  K  and  k  are  constant  and  Iq  is  set  equal 
to  zero.  In  ice,  K  and  k  can  vary  due  to  the  effects  of  brine 
pockets  as  noted  above  [Untersteiner,  1961]. 

Although  the  Maykut  and  Untersteiner  model  provides 
a  thorough  treatment  of  the  snow-ice  thermodynamic  sys¬ 
tem,  the  complexity  of  the  model  can  make  it  impractical  for 
modeling  large  areas.  Semtner  [1976]  simplified  the  Maykut 
and  Untersteiner  model  by  fixing  the  snow  and  ice  conduc¬ 
tivities.  The  salinity  piofile,  required  by  the  Maykut  and 
Untersteiner  model,  does  not  have  to  be  specified  in 
Semtner’s  version.  Internal  melting  is  attributed  to  an 
amount  of  penetrating  radiation  stored  in  a  heat  reservoir 
without  causing  ablation.  Energy  from  this  reservoir  keeps 
the  temperature  near  the  ice  surface  from  dropping  below 
freezing  in  the  fall.  Semtner’s  model  was  further  simplified 
by  assuming  linear  equilibrium  temperature  gradients  in 
the  snow  and  ice.  The  heat  flux  is  uniform  in  both  the  ice  and 
snow.  No  temperature  flux  exists  at  the  snow-ice  surface. 
The  surface  heat  balance  is  set  equal  to  zero  and  is  used  to 
solve  for  the  surface  temperature. 

The  third  category,  the  dynamic-thermodynamic  model, 
integrates  ice  motion  and  ice  groivth  and  decay  effects  into 
one  model.  Parkinson  and  Washington  [1979]  designed  one 
of  the  first  three-dimensional  dynamic-thermodynamic 
models.  Within  each  grid  cell,  the  model  was  broken  down 
vertically  into  four  layers:  a  mixed  layer  ocean,  an  ice  layer, 
a  snow  layer,  and  an  atmospheric  boundary  layer.  At  the 
ice-ocean  interface,  geostrophic  ocean  currents  were  used 
to  supply  ocean  stress  to  the  ice.  A  temporarily  invaniant 
dynamic  topography  and  a  constant  oceanic  heat  flux  were 
prosci  ibed.  Atmospheric  heat  fluxes  and  wind  stress  were 
applied  at  the  air-ice  interface.  The  ice  drift,  derived  from 
a  free-drift  formulation,  was  iteratively  corrected  to  insure 
the  existence  of  a  minimum  fraction  of  leads.  This  iteration 
does  not  conserve  momentum  and  can  damp  the  ice  drift. 

A  more  realistic  treatment  of  ice  dynamics  was  used  by 
Hibler  [1979],  This  model  uses  a  viscous-plastic  constitu¬ 
tive  law.  A  plastic  rheology  is  used  to  relate  the  strength  of 
ice  interaction  to  a  two-level  ice  thickness  distribution.  The 
ice  strength  P  is  defined  by 

P  =  P*h  (6) 

where  A  is  the  ice  compactness  or  concentration  and  P*  and 
C  are  constants.  Ice  concentration  is  defined  in  these 
models  as  the  fraction  of  a  grid  cell  that  is  covered  by  thick 
ice.  The  area  covered  by  thin  ice  or  open  water  is  1  —  A. 
Equation  (6)  allows  the  ice  strength  to  be  greater  in  regions 
of  ice  inflow  and  weaker  in  regions  of  ice  outflow. 


dA  d(uA)  d(i’A)  ,  n  c 

. — .  i  +  +(Diffu.s!cni  ibi 

dt  dx  dy 

where  S;,  and  are  source  and  sink  terms  accounting  for 
the  growth  and  melt  of  ice.  In  later  papers  le.g.,  Hibler, 
1980a],  these  growth  rates  were  calculated,  as  opposed  to 
prescribed,  from  the  heat  budget  balance  similar  to  that 
used  by  Parkinson  and  Washington  [19791 

Snow  cover  was  not  included  explicitly  in  the  Hibler  ice 
model.  Instead,  based  on  the  work  of  Bryan  et  ai.  1 1 975 !  and 
Manabe  etal. [1979], the  effects  of  snow  cover  we  re  a  pproxi  - 
mated  by  setting  the  ice  surface  albedo  to  that  of  snow  when 
the  surface  temperature  was  below  freezing  and  to  that  of 
snow-free  ice  when  the  surface  was  at  the  melting  point. 
Walsh  et  al.  [1985]  were  the  first  to  include  snow  explicitly 
in  a  long-term  simulation  using  daily  winds  and  air  tem¬ 
perature  data  from  the  National  Center  fo*-  Atmospheric 
Research  (NCAR)  for  the  years  1951  through  1980.  Inthis 
model,  both  snow  and  ice  were  composed  of  seven  thickness 
levels.  The  seven  levels  were  defined  by  taking  twice  the 
average  ice  or  snow  thickness  within  a  grid  cell  and  dividing 
itinto  seven  equal  levels.  The  growth  and  melt  of  the  ice  and 
snow,  based  on  the  heat  budget  balance,  was  calculated  at 
each  level  and  then  averaged  be  'k  into  one  thick  ice  growth 
rate.  Snowfall  rates  consisted  of  monthly  varying  climato¬ 
logical  rates  from  Maykut  and  Untersteiner  1 1971 J.  In  this 
system,  heat  from  the  heat  budget  balance  is  used  to  melt 
all  the  snow  before  it  is  used  to  melt  ice.  If  a  snow  cover 
exists,  the  albedo  is  set  to  0.80.  When  there  is  no  snow,  the 
albedo  is  0.65.  A  version  of  the  Hibler  ice  model,  similar  to 
that  used  by  Walsh  et  al.  [1985],  is  presently  used  as  the 
basis  for  the  Polar  Ice  Prediction  System  fPIPS),  the  U.S. 
Navy’s  numerical  sea-ice  forecast  model. 

22.2.2  Ice-Ocean  Models 

The  development  of  models  for  ice  covered  regions  was 
further  expanded  by  including  the  temporal  and  spatial 
variability  of  the  ocean.  This  was  accomplished  by  fully 
coupling  ice  to  ocean  models.  The  first  fully  coupled ,  three- 
dimensional  model  was  designed  by  Hibler  and  Biy  an 
[1984, 1987].  The  model  coupled  the  Hibler  ice  model  to  the 
Bryan-Cox  multilevel  ocean  model  [Bryan,  1969]  both 
dynamically  and  thermodynamically.  The  ocean  model  wa.-^ 
initialized  from  climatological  temperature.^i  and  .salinities 
and  weakly  constrained  to  those  values  with  a  three-year 
relaxation  period.  Tests  performed  using  this  model  high¬ 
lighted  the  importance  of  the  deep-ocean  heat  flux  into  the 
mixed  layer.  This  heat  is  great  enough  to  keep  parts  of  the 
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Barents  and  Greenland  Seas  ice-free  all  winter.  Semtner 
[1987]  used  a  similar  coupled  ice-ocean  model,  but  removed 
the  constraint  on  the  ocean  temperature  and  salinity.  He 
also  used  simplified  ice  dynamics  that  contained  only  bulk 
viscosity  (no  shear).  Semtner’s  results  verified  the  impor¬ 
tance  of  the  oceanic  heat  flux  on  keeping  the  marginal  ice 
zone  ice-free  in  winter.  However,  his  ice  thickness  distribu¬ 
tion  contained  thinner  ice  than  found  by  Hibler  and  Brj'an. 
Fleming  and  Semtner  [1991  ]  reduced  the  ice  strength  used 
by  Semtner  [1987]  and  obtained  a  much  more  realistic  ice 
thickness  distribution.  Riediir.jcr  and  Preiler  (1991  j  u.sed 
a  model  similar  to  the  Hibler  and  Bryan  model,  but  used 
daily  varying  forcing  from  the  Navy’s  operational  global 
atmospheric  model.  This  study  showed  that  the  diagnostic 
ocean  model  could  provide  a  useful  tool  for  forecasting. 

Additional  ice— ocean  models  have  been  developed  that 
show  the  importance  of  including  a  mixed  layer  in  the  ocean. 
One-dimensional  models  [Pollard  et  al.,  1983;  Lemke  and 
Manley,  1984;  Ikeda,  1985;  Lemke,  1987;  Bjork,  1989; 
Mellor  and  Kantha,  1989;  Riedlinger  and  Warn-Varnas, 
1 990],  two-dimensional  models  [Kantha  and  Mellor,  1989], 
and  three-dimensional  models  [Piacsek  ot  al .,  1 99’ }  have  all 
shown  distinct  changes  to  the  heat  and  salt  exchange  (and 
therefore  the  ice  growth  and  decay  rates)  when  a  mixed 
layer  is  included. 

These  dynamic-thermodynamic  ice  and  ice-ocean  mod¬ 
els  have  mainly  been  applied  to  large-scale  simulations 
using  grid  resolutions  on  the  order  of  1 00  km.  The  Hibler  ice 
model  has  been  applied,  however,  on  regional  scales  to  the 
Greenland  Sea  [Tucker,  1 982;  Preiler  et  al.,  1 990]  and  to  the 
Barents  Sea  (Preiler  et  al..  1989]  at  grid  resolutions  of 
approximately  20  km. 

Higher  resolution  models  (of  a  few  kilometers)  have 
often  been  used  to  look  at  processes  normally  associated 
with  the  region  near  the  ice  edge.  Roed  and  O’Brien  [1983] 
and  Hakkinen  [1986]  used  one-  and  two-dimensional  ice- 
ocean  models,  respectively,  to  study  upwelling  at  the  ice 
edge.  Ikeda  [1988]  used  a  three-dimensional  ice-ocean 
model  that  included  thermodynamics  to  investigate  upwell¬ 
ing  at  the  ice  edge.  Tang  [1991]  used  a  two-dimensional 
thermodynamic  ice-ocean  model  to  study  the  advance  and 
retreat  of  the  ice  edge  near  Newfoundland.  Smith  et  al. 
[1988]  used  a  three-dimensional,  two-layer  ocean  model 
coupled  to  a  Hibler-type  ice  model  to  investigate  the  behav¬ 
ior  of  an  isolated  ocean  eddy  within  the  marginal  ice  zone. 
Smith  and  Bird  [1991]  used  the  same  model  to  study  the 
interaction  of  an  ocean  eddy  with  a  jet  at  the  ice  edge. 

22.2.3  Ice-Atmosphere  Models 

Because  ice-atmosphere  interaction  cannot  be  divorced 
from  the  ocean’s  influence,  there  have  been  relatively  few 
studies  utilizing  interactive  models  of  only  the  atmo.spherc 
and  sea  ice.  The  studies  that  fall  into  this  category  have 
generally  been  one-  or  two-dimen.sional  experiments  with 
energy  balance  models.  Such  .studies  have  tended  to  focus 
on  the  effects  of  specific  proces.se.s.  For  example,  Lcdley 


[1988,  1991]  has  u.st-d  a  coupled  energy-balance  tea-ice 
model  to  examine  the  effects  of  presenbed  It  ad  frartiiin.s. 
siiowfiill,  and  a  highly  parameterized  iraiLSpori  of  sea  ice 
Ledley’s  model  i.s  one-dimensional  i. latitudinal;,  but  U  doe.c 
mcUide  a  partitioning  of  each  latitudinal  .'.one  into  land  and 
ocean.  Harvey  [19881  developed  a  more  cotniirehensive  sen 
ice  model  for  use  in  zonally  avera.ged  energy-lialance  cli¬ 
mate  siinu!ation.s.  Harveys  model,  which  includt.'.s 
parameterization.^  of  processes  such  a.s  surface  and  lateral 
melting,  lead.s,  advectioti,  and  snow  and  ice  ihickne.s.-i  ditt  n 
butions,  has  been  u.sed  to  study  variou.s  .sea  ict.*  feedback.s 
(albedo,  leads,  etc.).  Although  model-s  .such  a.s  these  are  so 
highly  idealized  th,al  their  relevance  to  the  actual  dimaU- 
system  is  unclear,  they  do  identify  high-leverage  porans 
eters  and  proce.sses  that  nicnl  further  .study  with  more 
sophisticated  models. 

Ice-atmosphere  models  have  also  been  used  to  study 
atmospheric  boundary  layer  proce.sse.s  and  interaction.-  at 
the  air-ice  interface.  Overland  [1 98.5  j.  for  example,  u.sed  an 
atmospheric  boundary  layer  model  to  addre.ss  momentum 
exchange  and  surface  drag  parameienzation.s  at  the  air-ice 
interface.  Bennett  and  Hunkins  1 1 986]  u.sed  a  two-dimen¬ 
sional  model  to  examine  atmospheric  boundary-layer  niodi- 
ficalions  during  advcction  over  an  inhomogencou.s  sea  ice 
cover.  Pease  and  Overland  II 9S4J  studied  the  dnfi  and 
mass  balance  of  sea  ice  in  the  Bering  Sea  by  using  an 
atmospherically  driven  sea  ice  model. 

22.2.4  Linkage  to  General  Circulation  Models 

A  survey  of  current  versions  of  the  genera!  circulation 
models  (GCM’s)  used  to  simulate  the  global  climate  reveals 
a  substantial  gap  between  state-of-the-art  sea  ice  models 
and  the  treatment  of  sea  ice  in  GCM's.  Most  climate  models 
treat  sea  ice  as  a  motionless  slab  characterized  by  a  single 
thickness  in  each  cell.  The  simulated  sea  icc  simply  accreie.s 
or  melts  in  response  to  a  deficit  or  surplus  in  the  .surface 
energy  budget.  The  treatments  that  permit  sea  ice  motion 
do  so  with  relatively  crude  parametorizations.  For  ex¬ 
ample,  the  model  of  Manabe  and  Stouffer  i  1 988  j  retains  the 
formulation  of  Bryan  [1969],  whereby  icc  moves  in  the 
direction  of  the  surface  ocean  currents  until  the  ice  reackc.- 
a  threshold  thickness  (4  m),  after  which  ice  motion  cea.-e.- 
Few  models  include  leads  and  open  water  area.-  within  the 
pack  in  spite  of  the  importance  of  these  element-  to  the 
atmospheric  simulation,  as  shown  in  variou.s  .sensitivity 
studies (e.g..[Simmonds and  Budd.  1990]).  The  model.- that 
do  include  leads  generally  prescribe  the  lead  fraction.-,  as 
discussod  below. 

The  inclusion  of  icc  transport,  lead.s,  and  open  water  in 
interactive  models  is  es.sential  to  the  realistic  simulat  ;on  of 
sea  ice  and  climate  for  several  rea.-on.s.  Areas  oft  bin  ice  and 
open  water,  continually  created  by  deformation,  are  the 
sources  of  nearly  all  the  new  ice  growth,  salt  rejection,  and 
heat  exchange  with  the  atmosphere.  The  flux  of  risible 
heat  from  the  ocean  to  the  almo.-phere  during  winter  r-  one 
to  two  order.s  of  magnitude  greater  over  thin  icc  and  open 


water  than  over  perennial  ice.  Maykut  [1978]  has  shown 
that,  in  an  areally  averaged  sense,  a  large  fraction  of  the 
heat  transfer  from  the  polar  oceans  to  the  atmospheric 
boundary  layer  during  winter  takes  place  over  young  ice 
(thickness  <0.5  m).  The  summer  melt  of  ice  and  the  heat 
storage  in  the  ocean  are  accelerated  dramatically  by  the 
absorption  of  solar  radiation  in  leads  and  other  low-albedo 
areas  such  as  new  ice. 

With  regard  to  ic?  tran.sport,  the  climatological  mean 
pattern  of  ice  motion  across  the  Arctic  Basin  and  through 
the  Fram  Strait  results  in  a  net  salinity  flux  to  the  continen¬ 
tal  shelves  off  Siberia.  These  shelves  are  regions  of  net  ice 
growth.  The  Greenland,  Iceland,  and  Norwegian  Seas  are 
regions  of  net  ice  melt.  The  net  gc  in  of  salt  by  brine  rejection 
in  the  shelf  waters  of  the  Arctic  Ocean  is  thought  to  drive  the 
thermohaline  circulation  of  the  Arctic  Ocean  [Aagaard  et 
al.,  1985].  The  export  of  sea  ice  into  the  North  Atlantic 
subpolar  seas  is  a  major  source  of  fresh  water  for  this  region, 
in  which  deep  water  formation  occurs  intermittently  and  is 
highly  sensitive  to  the  salinity-determined  stratification  of 
the  ocean  surface  layers  [Aagaard  and  Carmack,  19891. 

Although  leads  are  not  determined  by  ice  dynamics  in 
current  versions  of  major  GCM’s,  prescribed  or  simply 
parameterized  leads  hav?  recently  been  included  in  several 
climate  simulations.  The  United  Kingdom  Meteorological 
Office  (UKMO)  and  Australian  model,  fcr  example,  have 
been  run  with  seasonally  varying  prescribed  lead  fractions. 
Sensitivity  experiments  performed  with  the  Australian 
model  indicate  that  prescribed  lead  fractions  of  50%  in  the 
Antarctic  during  winter  cause  the  simulated  circulation  to 
resemble  more  closely  the  case  of  an  ice-free  southern  ocean 
than  the  case  of  an  ice-covered  southern  ocean  [Simmonds 
and  Budd,  1990].  Even  when  relatively  modest  lead  frac¬ 
tions  (several  percent)  are  prescribed  in  accordance  with 
observations,  the  atmospheric  sea  level  pressure  changes 
by  5  mb  in  some  high-latitude  regions.  An  alternative 
strategy  has  been  followed  in  the  Goddard  Institute  for 
Space  Sciences  (GISS)  II  model  (Raymo  et  al.  ,1990],  which 
computes  a  fraction  of  open  water  (leads)  as  fr  =  O.l/Zj^g, 
where  is  the  ice  thickness  in  meters.  The  latter  is 
currently  prescribed  as  1  m  at  all  times,  resulting  in  a  lead 
fraction  of  10%;  the  use  of  computed  values  of  Zj^  will 
permit  more  thorough  tests  of  this  parameterization.  This 
parameterization  of  leads  is  mentioned  here  merely  to 
indicate  the  relatively  simplistic  nature  of  such 
parameterizations  instate-of-the-artglobal  climate  models. 

1 1  should  be  noted  that  the  ice-ocean  interactions  involv¬ 
ing  salinity  (freshwater)  fluxes  will  not  be  reproduced 
correctly  by  coupled  models  until  the  models  include  ice 
transport.  While  ice  transport  is  simulated  well  by  atmo¬ 
spherically  forced  ice  models  (e.g. ,  PIPS),  a  successful  Simu¬ 
la  tion  of  ice  transport  in  coupled  GCM’s  will  require  an 
adequate  simulationof  the  surface  wind  field  and,  hence,  of 
the  sea  level  pressure  field.  There  is  no  evidence  yet  that 
atmospheric  GCM’s  reproduce  the  Arctic  pressure  patterns 
that  will  produce  the  key  features  of  ice  motion  observed  in 
theBeaufort  Gyre,  transpolar  drift  stream,  and  East 
Greenland  drift. 


22.2.5  Process  Models 

In  addition  to  the  model  examples  given  in  Section,-, 
22.2.1  through  22.2.4,  a  variety  of  other  phy.sical  inodeLs 
exist  to  represent  the  individual  components  of  sea  ice 
models,  including  detailed  thermodynamic-s  of  hea'  tr.ins- 
port,  turbulent  flu.x  models  at  the  air~<xean  interface, 
atmospheric  boundary  layer  processe.s,  radiative  iiiodeLs  of 
longwave  and  shortwave  radiation,  and  .simplified  oceanic 
mixed  layer  models. 

As  noted  earlier,  heat  transport  through  the  .snowpack  i.s 
ty^jically  treated  in  tenn.s  of  conduction  only,  v  ilh  one 
average  thermal  conductivity  u.sed  for  the  entire  .snow 
column  (e.g..  [Parkinson  and  Washington,  1979;  .Semlner, 
1976;  Maykut,  1 982]).  Colbc-ck  [  1 991  j  dj.<cu.s.se.s.some  of  the 
approaches  used  to  addrcs.s  conduction  as  well  a.s  turbulent 
fluxes  through  homogeneou.s  and  layered  snow-packs.  Pen¬ 
etration  of  the  snowpack  at  optical  wavelerunhs  i.s  a.ssurned 
tobenegligible[ParkinsonandWashington,1979;  '  .reated 
with  a  s.ngie  transmittance  value  <e.g.,  [Gabi.^on.  1987)  auo 
Section  22.2.1),  More  detailed  treatments  of  thenna!  and 
radiative  propertie.sof  the  snowpack  are  de.scnbed  by  Weller 
and  Schwerdtfeger  [1977]  and  Schwerdtfeger  and  Weller 
[1977],  Turbulent  fluxes  at  the  snow/ice— air  .surface  and 
open-water-air  surface  are  treated  simpH.sticaily  in  three- 
dimensional  ice  models,  with  a  single  energy  tran.sfer  coef¬ 
ficient  for  a  variety  of  different  surface  types.  Andreas  and 
Murphy  [1986],  Andreas  [1987],  Morris  !1989|,  and  Smith 
et  al.  [1990]  discuss  more  detailed  treatments  of  turbulent 
transferoverice  surfaces  and  leads.  Specific  processes  ofice 
growth  and  oceanic  modifications  in  refreezing  leads  are 
discussed  by  Bauer  and  Martin  [1983!  and  Kozo  [1983], 

In  addition  to  the  treatments  of  atmosphere— ocean-ice 
coupling  presented  earlier,  simplified  models  exist  with  the 
potential  to  improve  the  performance  of  the  existing  ice 
models.  Overland  et  al.  [1983],  Chu  [1987],  and  Overland 
[1988]  describe  models  of  the  atmospheric  boundary  layer 
over  the  ice  pack.  A  simplified  coupling  of  the  ice  cover  to  the 
boundary'  layer  is  described  by  Koch  1 1 98SJ.  Simple  repre¬ 
sentations  of  the  upper  ocean  applicable  to  ice  modeling  are 
described  by  Fichefet  and  Caspar  [1988]  and  Wood  and 
Mysak  [1989],  A  variety  of  parameterizations  to  improve 
the  accuracy  of  the  longwave  and  shortwave  fluxe.s  that  are 
input  to  the  ice  models  are  available,  ranging  from  simpli¬ 
fied  approaches  [Harvey,  1988;  Shine  and  Henderson- 
Sellers.  1985;  Ebert  and  CuiTy,  19901  to  more  complt-lc 
transmittance  models  [Kjteizys  et  al.,  1  988) 

22.2.6  ModelParameteruatto7isofFiclds.-\La}lah!cFro}7i 
Sotcllite  Data 

The  previous  section  contained  a  .summary  ol  'he  type.« 
of  numerical  models  that  presently  exi.-i  to  de.-enbe  ihv 
interaction  of  atmosphere,  ice,  and  ocean  I  n  order  to  better 
clarify  the  relationship  between  these  numerical  model.-; 
and  remotely  .sensed  data,  the  following  description  of 
existing  parameterizations  ofkey  model  fields  is  presenud 
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In  particular,  the  fields  that  are  or  will  soon  be  available 
operationally  from  satellite  data  are  ice  extent,  concentra¬ 
tion,  multiyear  fraction,  and  velocity. 

Ice  extent  is  determined  in  a  relatively  straightforward 
manner  from  models  containing  sea  ice  thermodynamics 
and/or  dynamics.  When  the  temperature  of  the  mixed  layer 
(or  upper  layer  of  the  ocean)  falb  below  freezing  (approxi¬ 
mately  —1 .8'^C)  in  a  particular  grid  cell,  sea  ice  forms  in  that 
grid  cell.  The  ice  is  subsequently  advected  into  other  grid 
cells  or  it  remains  in  its  area  of  formation;  melt  may  occur 
at  any  time.  Ice  extent  can  be  defined  as  the  equatorward 
limit  of  sea  ice  (of  any  concentration).  Alternatively,  ice 
extent  is  often  expressed  quantitatively  as  the  ocean  area 
poleward  of  this  limit  or,  e.g.,  in  the  NASA  sea  ice  atlases 
[Zwally  et  ai.,  1983;  Parkinson  et  al.,  1 987]  as  the  area  with 
ice  of  concentration  >15%.  The  areal  summation  is  per¬ 
formed  similarly  whether  one  is  evaluating  ice  extent  from 
obseiwational  (satellite)  data  or  from  model  output.  Slight 
discrepancies  may  be  introduced  by  the  spatial  resolutionof 
the  data  or  model  output  from  which  ice  ey*ent  is  evaluated. 

Ice  concentration,  which  provides  direct  information  on 
the  open  water  fractions  within  an  area  of  sea  ice,  can 
change  by  several  mechanisms  in  model  simulations:  ad- 
vection,  divergence  or  convergence,  and  freezing  or  melt. 
Small  changes  of  ice  concentration  can  also  be  caused  by 
diffusion  terms,  which  are  included  in  some  models  in  order 
to  mimic  the  effects  of  subgrid-scale  motions  or  to  provide 
numerical  stability.  Melt  can  affect  ice  concentrations  by 
changing  the  thickness  (to  zero)  or  by  melting  the  ice 
laterally  (e.g.,  in  leads).  In  models  containing  only  a  single 
ice  thickness,  melting  at  the  surface  can  change  the  concen¬ 
tration  without  melting  all  the  ice  in  a  grid  cell  if  the  single 
thickness  is  assumed  to  represent  the  mean  of  a  range  of 
thicknesses,  e.g.,  a  linear  distribution  of  thicknesses  be¬ 
tween  zero  and  twice  the  mean  (e.g.,  (Hibler,  1979]).  This 
strategy  has  not  been  adopted  in  GCM’s,  but  it  has  been 
used  in  various  stand-alone  ice  models  [Preller  and  Posey, 
1989;  Walsh  and  Zwally,  1990].  Changes  of  ice  concentra¬ 
tion  by  freezing,  on  the  other  hand,  are  typically  assumed  to 
create  ice  of  a  uniform  thickness  over  all  open  water  in  a  grid 
cell  unless  the  model  formulation  includes  a  prescribed 
minimum  fraction  of  leads  or  open  water. 

The  formulation  of  the  multiyear  ice  fraction  is  based 
partially  on  the  definition  of  multiyear  ice:  sea  ice  that  has 
survived  at  least  one  summer’s  melt  season.  Consequently, 
models  that  carry  multiyear  ice  as  a  variable  impose  an 
abrupt  transition  from  first-year  ice  to  multiyear  ice.  This 
transition  can  occur  either  at  the  time  when  the  mixed  layer 
temperature  drops  to  freezing  or  at  a  prescribed  date  (e.g., 
September  1  in  the  Arctic),  When  multiyear  ice  is  present, 
its  concentration  changes  by  the  same  processes  thatchange 
total  ice  concentration:  advection,  convergence  or  diver¬ 
gence,  diffusion,  and  melting,  (Multiyear  ice,  by  definition, 
cannot  form  instantaneously  by  freezing.)  Surface  heat 
budgets  are  computed  separately  for  the  first-year  and 
multiyearice  within  agridcell.  The  ice  rheology  determines 
the  extent  to  which  convergence  deforms  first-year  and 


multiyear  ice.  According  to  present  algorithms  1  Walsh  and 
Zwally,  ]  990],  convergence  deforms  solely  first-year  ice  as 
long  as  the  concentration  of  first-year  ice  exceeds  5% .  The 
fcnruilation  of  the  ice  strength,  which  ha.s  been  made  a 
rather  ad  hoc  function  of  the  ainountt;  of  flr.sl  year  and 
multiyear  ice,  needs  further  attention  in  modehs  that  distin¬ 
guish  several  types  of  ice.  Hiblcr’s  [1980al  multiple-thick¬ 
ness  formulation  provides  a  possible  point  of  departure  for 
formulations  of  the  strength  of  first-yea  rand  multiyear  ice 
mixtures. 

22.3  Model  Requirements  oe  SATt:LLnr:  Data 

Ice-ocean  models  require  obser\’ational  data  for  their 
initialization  and  validation.  Moreot'er,  the  use  of  observa¬ 
tional  data  for  the  forcing  of  ice-ocean  models  provides  an 
objective  means  for  model  testingand  assessment  when  the 
models  are  run  in  a  decoupled  (from  the  atmosphere )  mode. 
The  general  lack  of  in-situ  measurements  from  ice-co  vercd 
waters  makes  satellite  remote  sensing  the  cornerstone  of 
data  assimilation  schemes  required  for  ice  modeling. 

The  manner  in  which  data  are  used  for  model  initializa¬ 
tion  depends  on  which  of  three  types  of  ice  model  s  i  m  ul  aiion 
is  being  performed:  short  term  (days),  long  term  wath 
periodic  forcing,  and  long  term  with  aperiodic  forcing.  An 
example  of  short-term  simulation  is  the  use  of  an  opera¬ 
tional  model  to  forecast  sea  ice  at  ranges  of  several  days 
(e.g.,  PEPS).  For  such  simulations,  an  accurate  initializa¬ 
tion  corresponding  to  the  actual  sea  ice  state  is  cruci  a!  to  the 
success  of  the  forecast.  The  second  type  of  model  integra¬ 
tion,  a  long-term  periodic  simulation,  spans  at  least  several 
years  with  an  annual  cycle  of  forcing  that  repeats  itself 
exactly.  In  such  cases,  the  ice  model  should  achieve  equilib¬ 
rium  if  run  for  a  sufficiently  long  time.  This  type  of 
simulation  is  useful  for  diagnostic  studies  of  the  ice— ocean 
system  (e.g.,  [Hibler  and  Bryan,  1987;  Semtner,  1987];  and 
for  assuring  that  the  ice  model  is  free  of  long-term  biases  or 
model  drift  prior  to  interactive  coupling  isuch  as  to  an 
atmospheric  (jCM).  Provided  that  sufficient  computational 
resources  are  available,  a  realistic  initialization  is  not 
essential  for  such  simulations  because  the  model  should 
eventually  achieve  equilibrium  with  its  forcing.  The  third 
type  of  simulation,  the  long-term  aperiodic  run,  is  used  to 
obtain  an  interannually  varying  history  tape  of  the  .-^ea  ice 
state.  Examples  include  the  hindcast  simulation.s  of 
Rothrock  and  Thomas  [1990],  W’alsh  and  Zwally  1 1 990],  and 
Riedlinger  and  Preller  [1991],  As  in  the  periodic  simula¬ 
tions  described  above,  the  mode!  can  be  spun  up  to  the 
initial  state  of  the  desired  date  by  integrating  ( periodically 
or  aperiodically)  through  several  annual  cycles  I^or  this 
reason  and  because  the  long-term  statistics  of  the  simula¬ 
tion  are  of  primary  interest,  a  precise  initial  .state  is  not 
essential  to thesuccessofthesimulation,  Neverthele.s.s,  the 
initial  state  must  conform  reasonably  closely  to  the  corre¬ 
sponding  actual  state  if  the  first  several  year.s  are  not  to 
contaminate  the  verification  statistics.  Thus,  the  spi  n-u  p  ( if 
such  a  simulation  must  be  achieved  in  such  a  way  that  it 


ensures  temporal  homogeneity  in  the  history  tape’s  corre¬ 
spondence  to  observations.  One  may  argue  that  the  initial¬ 
ization  strategy'  involves  more  subtleties  in  this  type  of 
simulation  than  in  the  other  two. 

The  initialization  and  validation  of  ice  models  require 
observ'ational  data  on  sev'eral  key  sea  ice  variables:  concen¬ 
tration,  thickness  and  ice  type,  drift  velocity  and  associated 
kinematic  quantities,  and  surface  properties  (albedo,  snow 
cover,  temperature, etc.).  Whileother propertiessuchasice 
salinity  and  temperature  profiles  may  be  important  for 
diagnostic  studies  with  ice-ocean  models,  we  will  focus 
primarily  on  the  four  variables  listed  above.  These  proper¬ 
ties  are  also  variables  for  which  satellites  provide  useful 
information. 

22.3.1  Key  Variables  Used  for  Model  Initialization  and 
Validation 

22.3.1.1.  Ice  concentration.  Ice  concentration  data  im¬ 
plicitly  include  information  on  several  quantities  that  are  of 
central  importance  to  studies  of  the  large-scale  distribution 
and  variability  of  sea  ice;  fractional  coverage  of  an  area 
(ideallj'  within  different  ice  thickness  ranges;,  ice  extent, 
and  the  fraction  of  open  water  (e.g.,  leads)  within  the  pack. 
For  forecasting  purposes,  ice  extent  may  be  the  primary 
concern.  For  coupled  model  studies  of  the  Arctic  or  Antarc¬ 
tic  surface  energy  balance,  the  lead  fraction  may  be  a 
primarj’  concern  because  the  small  areas  of  leads  and  open 
water  effectively  govern  the  exchanges  of  heat  and  moisture 
between  a  pack-ice  surface  and  the  atmosphere.  Different 
sensors  provide  the  optimal  information  on  ice  concentra¬ 
tion,  extent,  and  lead  fractions. 

The  most  useful  and  homogeneous  source  of  ice  concen¬ 
tration  data  for  large  (50  km  x  50  km)  areas  is  the  passive 
microwave  data  from  polar  orbiting  satellites:  the  Electri¬ 
cally  Scanning  Microwave  Radiometer  (ESMR),  the  Special 
Sensor  Microwave  Imager  (SSM/I),  and  the  Scanning  Mul¬ 
tichannel  Microwave  Radiometer  (SMMR).  Since  1978, 
these  satellites  have  provided  nearly  continuous  coverage 
at  1-  to  2-day  intervals  and  50-km  resolution  with  more 
scattered  coverage  being  available  since  1972.  The  recent 
availability  ofgridded  passive  microwave  dataon  CD-ROM 
(through  the  National  Snow  and  Ice  Data  Center)  has 
greatly  facilitated  the  use  of  derived  ice  concentrations  by 
the  sea  ice  community. 

The  accuracy  of  any  data  used  for  model  verification  or 
data  assimilation  is  a  key  consideration  in  determining  its 
usefulness.  Errors  in  the  total  ice  concentration  (first-year 
plus  multiyear  ice)  and  ice  edge  location  as  derived  from 
passive  microwave  data  are  fully  discussed  in  Chapters  4, 
10,  and  11.  The  estimated  errors  for  total  ice  concentration 
generally  range  from  5%  [Cavalieri  et  al.,  1984;  Swift  and 
Cavalieri,  1985;  Steffen  and  Maslanik,  1988]  to  about  15% 
[Steffen  and  Schweiger,1990)  using  general  versionsof  the 
NASA  Team  algorithm,  with  errors  as  low  as  3%  when  the 
algorithm  is  tuned  for  local  conditions.  Absolute  accuracies 
remain  uncertain  and  somewhat  controver.sia!  due  to  the 


difficulties  of  comparing  different  data  lype>.  and  difleri,  ni 
scales.  Total  icc  concentration  estimale.s  are  relatively 
unaffected  by  surface  melt  on  a  regional  ba.si.s,  except  when 
melt  ponds  are  present  (e.g.,  [Steffen  and  Maslanik,  1 986]  t 
Gloerscn  and  Campbell  11 988]  describe  some  more  specific 
effects  ofsurface  melton  retrievals  of  ice  concent  ration  and 
ice  type.  In  addition  to  melt  effect.s,  sources  of  error  includir 
the  effect  ofsurface  wind  on  open  water,  atmospheric  water 
content,  changes  in  the  ice  surface  (snowpack,  Hooding, 
etc.),  and  the  pre.sence  of  radiomctrically  thin  ice  not  ac¬ 
counted  for  by  the  algorithm. 

For  broad  depictions  of  the  total  ice  concentration  fit -Id.s 
in  large-scale  model.s  (with  resolution-s  on  the  order  of  100 
km),  the  accuracy  and  resolution  of  the  pa.s.sive  microwave 
data  are  certainly  adequate,  except  perhaps  during  late 
summer  when  substantial  melt  ponding  is  pre.seiU.  How¬ 
ever,  mesoscale  modcLs  may  require  finer  resolution  and 
greater  accuracy  (for  the  determination  of  the  kKation  of 
leads  and  polynas  in  the  pack)  for  initialization  and  valida¬ 
tion  at  the  10-km ‘scale. 

The  concentrat  ions  derived  from  satellite  passive  micro- 
wave  data  provide  part  of  the  input  to  the  Navy  and 
National  Oceanic  and  Atmospheric  .Administration  (NO.A.^) 
Joint  Ice  Center’s  (JIC)  weekly  analyses,  which  have  been 
used  in  the  initialization  of  the  PIPS  model  and  in  the 
compilation  of  verification  statistics  from  long-term  simu¬ 
lations  [Walsh  et  al.,  1985],  The  JIC  depictions  of  the  ice 
edge  have  also  utilized  information  from  the  radar  altim¬ 
eter  on  the  Geosat  satellite.  Because  the  shape  of  a  radar 
altimeter’s  return  pulse  is  significantly  different  over  sea 
ice  and  open  water,  the  altimeter  permits  the  resolution  of 
the  ice  edge  to  within  10  km,  which  is  finer  than  the 
resolution  of  the  digitized  JIC  ice  charts  and  considerably 
finer  than  the  resolution  of  the  passive  microwave  concen¬ 
tration  grids.  For  simulations  with  me.soscale  models, 
radar  altimeter  data  can  serve  a  useful  purpose  in  the 
initialization  and  validation  of  the  ice  edge. 

Fine  resolution  imagery  is  needed  to  resolve  the  indi¬ 
vidual  small  areas  ofopen  water  Heads  and  polynas)  within 
the  pack  ice.  Together  with  areas  of  thin  ice.  these  open 
water  areas  are  the  primary' dete rminantsof  the  icestrength. 
the  ocean-atmosphere  heat  exchange,  and  the  evaporation 
from  the  ocean.  Because  the  open  water  areas  often  have 
scales  of  tens  to  hundreds  of  meters,  individual  open  water 
areas  cannot  be  resolved  by  current  passive  microwave 
satellite  imagery,  although  such  subresolution  feature.?  are 
included  in  the  total  open  water  area  within  the  pa.ssive 
microwave  field  of  view.  Moreover,  the  uncertainties  of 
several  percent  in  the  passive-microwave-derived  concen¬ 
trations  are  comparable  to  the  variations  in  the  .-mall 
percentages  of  thin  icc  and  open  water  that  typically  occur 
in  pack  ice(e.g.,  [Maykut,  1982];.  Both  visible  and  infrared 
imagery  can  provide  fine-ro.sotutiond  0  meter.s  to  I  kilome¬ 
ter)  observations  of  ice  concentration,  ice  extent,  lead  loca¬ 
tions,  and  floe  characteristics.  Vi.sible  and  reflected  infra¬ 
red  imagery  are  limited,  however,  by  both  cloud.i  and 
darkne.s.s  while  thermal  infrared  imageiy  i.s  limited  by 
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cloud  cover.  The  most  promising  source  of  all-weather,  fine- 
resolution  data  is  the  satellite-borne  Synthetic  Aperture 
Radar  such  as  those  on  board  ERS-1  (from  vvliich  data  are 
currently  becoming  available),  JERS-1,  and  Radarsat. 
Nominal  resolution  from  the  JERS-1  SAR,  for  example,  will 
be  30  ra. 

The  accuracy  of  ice  concentration  from  SAR  data  is 
discussed  in  Chapter  6.  Generally,  errors  in  SAR-derived 
ice  concentration  occur  primarily  when  thin  ice  is  present, 
since  thin  ice  can  have  backscatter  characteristics  similar 
to  that  of  open  water.  However,  a  comparison  of  ice  concen¬ 
trations  calculated  from  visible-band  aerial  photographs, 
passive  microwave  acquired  from  aircraft,  and  SAR  gener¬ 
ally  agreed  to  within  15'7r,  with  ice— water  contrast  and  ice 
signature  variability  accounting  for  the  error  [Burns  et  ah, 

1987] . 

Considerable  potential  exists  for  the  use  of  SAR  data  in 
model  parameterization  studies  directed  at  the  realistic 
simulation  of  the  thin  ice  and  open  water  areas  within  pack 
ice.  Mesoscale  models,  in  particular,  will  benefit  from  the 
order  of  magnitude  improvement  in  resolution  made  pos¬ 
sible  by  SAR.  In  addition  to  uses  by  the  ongoing  ice 
modeling  efforts  summarized  here,  such  information  may 
serve  as  a  stimulus  to  the  coupling  of  sea  ice  and  atmo¬ 
spheric  models  through  realistic  parameterizations  of  sur¬ 
face-atmosphere  heat  exchange. 

One  additional  use  of  satellite-derived  ice  concentrations 
is  the  specification  of  lead  fractions  in  climatic  models  as 
noted  in  Section  22.2.4.  In  the  Australian  global  climate 
model  [W.  Budd,  personal  communication,  1991]  for  ex¬ 
ample,  the  specified  lead  fractions  are  the  seasonally  vary¬ 
ing  fractions  of  open  water  derived  from  the  microwave 
radiances  of  the  ESMR  and  SMMR  instruments. 

22.3.1.2.  Ice  thickness  and  ice  type.  The  spatial  and 
temporal  variations  of  ice  thickness  are  notoriously  difficult 
to  quantify.  The  only  data  on  sea  ice  thickness  have  been 
obtained  from;  { 1 )  holes  drilled  through  the  ice,  a  logistically 
difficult  and  expensive  proceuure,  (2)  submanne  sonar 
data,  much  of  which  is  classified,  providing  snapshots  for 
single  year-months  along  specific  cruise  tracks;  and  (3) 
single-point  time  series  obtained  from  moored  acoustic 
sensors,  which  have  come  into  use  only  during  the  last  few 
years.  Sources  (1 )  and  (2)  have  permitted  the  computation 
of  a  seasonal  climatology  of  Arctic  ice  thickness  (eg.,  1  Bourk“ 
and  McLaren.  1992])  used  to  validate  ice  models  in  a 
general  way,  but  they  have  provided  little  information  on 
interannua!  variability. 

There  are  no  immediate  prospects  for  the  direct  measure¬ 
ment  of  ice  thickness  by  satellite,  although  relationships 
between  ice  thickness  and  passive  microwave  polarization 
ratios  have  been  demonstrated  (e.g. ,  [Steffen  and  Maslanik, 

1988] ).  The  estimates  of  multiyear  ice  concentration  ob¬ 
tained  from  passive  microwave  data  can  provide  crude  and 
quantitative  indicators  of  ice  thickness,  but  such  indicators 
are  inadequate  for  model  initialization  and  verification. 
Thus,  there  will  most  likely  be  a  continued  reliance  on  the 


fragmentary  thickne.s.s  measurement.';  obtained  from  !ii<- 
surface  and  subsurface  .sources. 

Several  models  c.xplicitlv  include  the  concentration  and 
tliickness  of  multiyear  ice  in  their  formuialion  IRothruck 
and  Thomas,  1 990;  Walsh  and  Zwally,  1  >90;.  The  re.sults  of 
the  latter  model  show  some  agreement  with  the  interanmial 
fields  of  multiyear  ice  concentration  derived  from  SM.MR 
How'ever,  uncertainties  in  the  SMMR-denved  coni  enira 
tions  of  multiyear  ice  are  substantial,  and  u.seful  informa¬ 
tion  on  multiyear  icc  cannot  be  oht.iined  during  the  .spring 
and  summer  months  due  to  changes  in  tlie  dielectric  prop¬ 
erties  of  the  ice  pack  with  the  onset  of  melt  i  see  Chapter  -1  o 
During  winter,  average  accuracy  for  ice  type  h.is  been 
estimated  to  be  within  about  20';  ICavalieri  et  a!..  Dd-i ;. 
although  the  error  may  be  considerably  larger  even  in 
winter  [Comiso,  19861.  For  these  reasons,  the  Ktilman  filler 
approach  of  Rothrock  and  Thomas  provides  an  atiractiv  e 
means  for  assimilating  the  limited  .satellite  data  into  a  nior* 
temporally  homogcneou.s  framework.  Tims  framework  of 
fers  the  potential  for  a  more  meaningful  validation  of  tlv 
time-varying  fields  of  multiyear  ice  coverage  land  ice  thick¬ 
ness)  .simulated  by  the  more  traditional  sea  ice  models, 

22.3.1.3.  Ice  velocity  and  associated  kinematic  variables. 
Because  ice  velocity  has  a  relatively  short  time  scale  i  da>  .s  o 
accurate  initialization  of  ice  motion  is  important  primarily 
for  short-term  simulations.  For  long-term  .simulations, 
however,  the  mean  patterns  of  ice  drift  determine  the 
transports  of  mass  and  salinity.  The  validation  of  the 
simulated  drift  is  therefore  an  important  element  of  mode! 
simulations  of  long  periods  (sea.sons  or  longer). 

The  primary  source  of  data  on  time- varying  ice  velocities 
has  been  the  collection  of  position  measurements  of  Arctic 
buoys  and  drift  stations.  The  buoy  data  have  been  objec¬ 
tively  analyzed  on  a  daily  ba.sis  si  iice  1979  by  the  University 
of  Washington’s  Polar  Science  Center.  These  velocilie.s 
have  been  used  to  verify  model  simulations  of  ice  velocity 
(e.g.,  (Serreze  etal.,  1989]).  The  .satellite  pas.sive  microwave 
delta,  although  way  coar.se  in  resolution  to  identify 
individual  ice  floes,  nonetheless  permit  the  derivation  of 
approximate,  generalized  velocity  vectors  on  the  basis  of 
changes  of  multiyear  ice  concentration  over  poriod.s  of 
several  months  to  a  season.  Such  derivations  arc  based  on 
the  supposition  that  multiyear  ice  behaves  somewhat  like  a 
tracer  in  subft-eezing  regions  poleward  of  the  ice  edge.  For 
example,  velocities  denvccl  from  ciianges  ui  iijtr  ui 

the  multiyear  ice  edge  have  been  found  to  be  spatially 
coherent,  interannually  variable,  and  consistent  with  iht> 
fields  derived  from  a  sea  ice  mode!  [Zwally  and  Walsh. 
1987]. 

The  most  promising  .source  of  information  on  individual 
floe  ice  velocities  is  SAR  imagery  as  discussed  in  Chapter 
1 8.  Pattern  recognition  techniques  can  bo  applied  to  .S,\R 
images  of  the  same  region  at  intervals  of  .several  days, 
thereby  providing  kinematic  vector  imagc.s  of  the  net  u  e 
motion  during  the  intervening  time  period  i  e.g  ,  |  Curlander 
etal.,  19851).  Emery  et  al.  11991  i  used  a  .similar  approach 


with  Advanced  V'ery  High  Resolution  HadionietoriAVHIiK) 
data  to  obtain  ice  velocity  vectors  for  the  F ram  Strait  region; 
these  vectors  were  conipared  directly  ■.■.■‘h  model-derived 
ice  velocities.  The  use  of  data  in  this  manner  is  likely  to 
provade  the  most  precise  means  of  initializing  and  validat¬ 
ing  ice  velocities  from  meso.'cale  models  and,  perhaps,  also 
of  validating  the  integrated  measure  of  the  divergence  and 
convergence  over  grid  cell  areas  of  larger  scale  models.  U 
should  be  noted  in  this  regard  that  the  time  difTerences  of 
multiyear  ice  concentrations  derived  from  sattdlitc  passive 
microwave  data  also  provide  potentially  useful  information 
on  sea  ice  convergence  and  divergt-nce  over  larger  areas 
(hundreds  of  kilometers  !,  as  discussed  in  Section  ‘i’i.-l.'i. 

22.3.2  Additional  Satellite  Data  for  Future  Use  tn  Mode's 

22.3.2.1.  Albedo  and  ice  surface  temperature.  One  of  the 
critical  thermodynamic  properties  in  a  .sea  ice  model  simu¬ 
lation  is  the  surface  albedo.  While  the  albedo's  broad 
dependencies  on  the  surface  state  (snow,  bare  ice.  and 
puddled  ice)  are  generally  known,  the  ;.paiial  and  temporal 
variations  ofsurface  albedo  in  ice-covered  regions  are  poorly 
documented.  In  view  of  the  high  leverage  of  the  surface 
albedo,  especially  in  sea  ice  simulations  that  addres.$  cli¬ 
mate  change,  the  parameterization  and  validation  of  sur¬ 
face  albedo  merit  increased  attention. 

Albedos  over  ice  cover,  on  scales  suitable  for  model  input, 
have  been  both  mapped  from  digital  AVURR  imagery  (c.g., 
[Rossowetal.,  1989]J  for  short  time  periods  and  interpreted 
from  Defense  Meteorological  Satellite  Program  (DMSPi 
Operational  Line-Scan  System  (OLS)  transparencies 
[Robinson  et  ah,  1986;  Scharfen  et  ah,  1987]  for  about  a  ten- 
year  period.  These  DMSP  albedos  were  tested  in  an  ice 
model  by  Ross  and  Walsh  [1987 1.  Robimson  et  ah’s  results 
cover  only  the  central  Arctic  during  the  summer  melt 
season  (May  through  August),  but  do  provide  quantitative 
depictions  of  the  interannual  variations  over  areas  of  1 000 
'km  X  1000  km.  Unfortunately,  the  variation  of  regional 
surface  albedo  is  intertwined  with  the  distribution  of  open 
W'ater  within  the  region.  Model  bia.ses  regarding  open 
water  fractions  may  overwhelm  the  effects  of  variation  in 
the  snow-ice  surface  albedo.  A  need  therefore  exists  for  a 
more  complete  utilization  of  high-resolution  satellite  imag¬ 
ery  to  address  the  validitj^  of  model  treatments  of  surface 
albedo  within  a  grid  cell.  It  is  likely  that  combined  usage  of 
SAR  and  high-resolution  visible  imagery  will  be  required 
for  a  meaningful  as.sessment  of  the  parametric  needs  rela¬ 
tive  to  surface  albedo. 

The  energy  balance  calculations  necessary  to  estimate  ice 
surface  temperature  for  ice  growth  and  decay  in  models  are 
substantial  .sources  of  uncertainty.  Direct  mea.surement  of 
ice  surface  temperatures  would  be  particularly  valuable  for 
ice  prediction  models.  Comiso  |1983i  used  data  from  the 
temperature-humidity  infrared  radiometer  (THIR)on  board 
Nimbu.s  7  to  estimate  ice  surface  temperatures  covering 
nine  months  in  the  Arctic.  Schlue.ssel  and  Grassl  [1990] 
describe  the  use  of  AVHRR  to  retrieve  water  and  ice  surface 


temperatures  in  .•Vntarclic  pol  vin  a.-.  Key  and  l.'arrv  !  _ 

describe  the  retric\  .-il  of  ice  surface  teinpi.-raturi,::;  and  al’iie 
dos  for  a  portion  of  the  Arctic  using  .AVHRR  imagi-t  N’  .An 
example  of  the  calculation  of  ice  alliedos  and  temper.)!  u.'-i' 
from  AVHRR  for  operalional  ice  foreca.siing  i.s  discus.-  cd  (>y 
Condal  and  Le  llOS-l  j.  With  the  e.Kception  of  the  ail»  do.- 
mapped  from  D.MSP  OL.S  data,  no  albedo  or  timperaiure 
data  set  has  been  prejiared  for  the  .Arctic  that  cover:-  more 
than  several  months.  Hender.son— Seller.-,  and  U  lisou  , !  ds,-; ; 
review  some  older  albedo  climalologie.s  .suitable  for  mouei 
input.  Work  ns  underway  to  refine  the  inelhodokigy  to  niaj,- 
ICC  surface  temperature.s  la.c  well  as  cloud  cover,  from 
.AVHRR  a.s  part  of  the  Eartii  Ob.servmg  Sy.stem  Polar 
Kxchange  at  the  Sea  Surface  i  program 

22.3.2.2.  Snoiv  cover.  Snosv  cover  Lscle.iri  v  a  kiyv 

ofthc  albedo  formulation  discussed  above.  Snow  dejitii  a!;'o 
alTects  the  wintertime  energy  budget  over  .sea  ice  Wsi)!) 
snowfall  isani.exlernaliy  pre.scnbed  quantity  in  theenrite;-:'. 
of  ice  modeling,  there  is  ncvertheli.'s.s  a  need  for  rno.'-e 
information  on  the  di.stribution  of  .snow  depth  within  ihe 
grid-cell  areas  of  models.  The  extent  to  which  altimeti  r.s 
and.''or  .S.AR  can  provide  such  data  is  unclear.  Atleinpl.s  i" 
use  passive  microwave  (SMMR)  data  to  estimate  ,snow 
depth  over  land  have,  however,  met  with  .some  succe.s.s  for 
specific  regions  [Chang  cl  al,  1987).  ,Adiscu.<,sionofihe  u.s<- 
of  passive  microwave  data  to  determine  snow  depths  over 
sea  ice  is  presented  in  Chapter  16.  Advances  in  studie.s  of 
snow  depth  in  regions  covered  by  sea  ice  are  needed  to 
permit  modelers  to  determine  whether  biases  in  .sinuiiated 
sea  ice  (c.g..  thickne.s.sj  can  bo  attributed  to  bia.se.s  in  the 
prescribed  snowfall  or  .snow  depth. 

22.3.2.3.  Sea  surface  temperature.  Sea  surface  tempera¬ 
tures  suitable  for  ice  modeling  have  been  availaliie  fur 
several  years  and  are  derived  from  ship  observation.s  'e.g., 
the  comprehem.sivo  ocean-atmosphere  data  set.  CO.-VDSi 
and  satellite  imagery  such  as  AVHRR.  including  the  opera¬ 
tional  generation  of  sea  surface  temperatures  by  .N'O.A-V 
NESDIS  'National  Environmental  Satellite.  Data  and  In¬ 
formation  Seiwice),  the  multichannel  sea  surface  tempera¬ 
ture  fields  'o.g..  [McClain  et  al..  1985;  Barton.  2.985]!  and 
from  sounder  data  (e  g,,  [Sus.skind  and  Reuter.  19S5i>. 
Minnell  [1990]  reviewed  .some  of  the  attributes  of  sea 
surface  temperatures  measured  from  space,  including  the 
problem  of  determining  skin  versus  bulk  temperature.-; 
(e.g.,  (Scliiue.ssei  ei,ul..  19S7 ; Sciducsscl  and  f  r 1 990 ■  t 
While  the  noi.se  level  of  individual  temperature  mea.sure- 
mcnt.s  from  AVHRR  is  about  0.1  K,  the  unmodeled  effect.-  of 
the  intervening  atmosphere  reduce  thi.s  accuracy  I  he 
Along- Track  Scanning  Radiometer,  carried  on  board  ERS- 
1,  i.s  designed  wdtli  both  a  lower  noise  level,  a.s  well  a.s  the 
ability  to  observe  two  path.s  through  the  almosphnn-.  to 
achieve  more  preci.se  and  accurate  .surface  tempcraturi's, 

22.3.2.4.  Cloud  cover.  Farameterization.s  of  radiative 
fluxes  in  ice  models  can  benefit  from  reliable  cloud  cover 
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observations.  Long-term  cloud  cover  statistics  derived  from 
satellite  data  for  the  polar  regions  include  the  three-dimen¬ 
sional  nephanalysis  carried  out  by  the  U  .S.  Air  Force  [Fye, 
1978;  Tian  and  Curry,  1989)  ♦he  International  Satellite 
Cloud  Climatology'  Prograri.  >  'iCCPj  [Kossow  el  al.,  1985; 
Rossow  and  SchifTer,  1  Sy  ’ and  analyses  ofTHIR  imagery 
IStowe  et  al.,  1988'  .enderson-Sellers  [1984)  describes 
some  of  these  climatologies,  as  well  as  cloud  climatologies, 
from  surface  observations.  Someofthe  problemsof compar¬ 
ing  clou  j  climatologies  based  on  satellite  data  versus  sur¬ 
face  observations  are  reviewed  by  Henderson-Sellers  and 
•'  cGuffie  11990],  and  dilTercnt  satellite-derived  cloud  sta¬ 
tistics  for  the  Arctic  are  discussed  by  McGuffie  et  al.  [  1 988], 
Serreze  and  Rehder  11990)  describe  a  climatology'  for  June 
spanning  1 2  years  of  DMSP  interpretations  for  the  western 
Arctic.  Barry  et  al.  [1987]  give  DMSP-denved  cloud  statis¬ 
tics  for  the  Arctic  for  April  through  June  1979  and  1980. 
Key  and  Barry  [  1 989]  and  Ebert  [1 989 1  describe  some  ofthe 
problems  of  automated  mappingof  polar  clouds  and  AVHRR- 
based  cloud  mapping  schemes  designed  specifically  for  the 
polar  regions. 

Outgoing  longwave  radiation  fluxes  are  available  di¬ 
rectly  from  AVHRR  data  [Gruber  and  Jacobowitz,  1985), 
High-Resolution  Infrared  Radiation  Sounder  data  [Ellingson 
et  al.,  1989],  and  the  broadband  channels  of  the  Earth 
Radiation  Budget  Experiments  (ERBE  )(e.g.,[Kopia,  1986; 
Jacobowitz  et  al.,  1984;  Barkstrom  and  Smith,  1986; 
Ramanathan  et  al.,  19891).  Relatively  long  time  series  of 
these  data  are  available  at  different  spatial  and  temporal 
scales  for  the  polar  regions. 

22.0.2.5.  Atmospheric  properties.  Atmospheric  proper¬ 
ties  (water  vapor,  liquid  water,  and  temperatures)  are 
determined  operationally  using  sounders  as  well  as  passive 
microwave  imagery.  Sounders  provide  information  on  the 
vertical  distribution  of  atmospheric  conditions,  while  pas¬ 
sive  microwave  images  provide  integrated  values  for  water 
content.  Calculations  are  complicated  over  ice  and  snow 
due  to  the  high  and  vairiable  emissivity  of  the  surface. 

Surface  air  temperatures  used  in  ice  model  simulations 
are  generally  prescribed,  either  from  objective  analyses  of 
station  data  or  from  atmospheric  model  output.  The  effects 
of  errors  and  biases  in  these  temperatures  need  assessment 
in  much  the  same  way  as  do  the  effects  of  errors  in  the 
albedo  and  snow  depth  as  discussed  above.  The  use  of 
.satellite  data  for  such  assessments  does  not  appear  to  have 
been  fully  exploited,  although  the  products  of  the  ERBE  and 
ISCCP  efforts  may  be  useful  in  this  regard. 

Finally,  wmd  speed  over  open  ocean  is  routinely  deter¬ 
mined  from  passive  microwave  imagery  and  can  be  mapped 
from  SAR  data.  Speed  and  direction  can  be  mca.sured  by 
scatterometers.  Over  sea  ice,  no  direct  way  of  measuring 
surface  winds  from  satellites  is  available,  although  cloud 
niotion.sLave  been  u.sed  to  map  upper-level  windsoverpolar 
regions  (Turner  and  Warren,  1989], 


22.4  Ex.-\.mi’les  ok  thf,  I.NTKRAn’io.N  ok  .Moohi.s  W::s!  !)  .:  > 
22.4.1  Model  Urilala  Comporisoti  Studies 

Comparisons  with  satellite  data  have  pr  ovided  one  oj ;  he 
major  vehicles  for  recent  assessments  of  sea-Ki-  model 
performance.  As di.scussed below,  some  model  rtsshi:.  ioi’.e 
provided  insights  to  the  validity  of  saleiiite-denvi.  d  ileid:-  >)! 
variables  that  are  only  indirectly  measurable 

Comparisons  between  model  output  and  .satcilue  liat.i 
generally  require  systematic  compilation.^  of  salelhtc  mea 
surements  of  one  or  more  of  the  vanable-s  .simulated  by  a 
model.  An  outstanding  example  of. such  a  compilation  i  - 1  fie 
record  ofthe  sea  ice  extent  and  concentration  derived  fn irn 
satellite  passive  microwave  mcaburcments(ESM.H..8M. 'JR. 
andSSM/I).  The  passive  microwave  data  a^eal.soa  prirna!-. 
input  to  the  JIC’s  weekly  ice  charts,  which  ht."e  been  used 
for  a  host  of  additional  modei-to-data  compari.son-  The 
availability  of  the  passive  microwave  grids  of  ice  coneen  t  ra 
tion  from  1973  onward  provides  a  climatology  a.~  w,  !i  as  a 
15- to  18-ycarsampleof  interannual  nuctualion.s  ai-.out  the 
seasonal  means. 

Ajnong  the  first  model  studie.s  to  have  utilized  data  on  ice 
extent  derived  from  passive  micrc'.vave  imagery  were  those 
ofParkinson  and  Henman  11980)  and  Lemko  el  al  il9s'H 
In  the  latter  study,  a  stochastic-dynamic  mode!  wa.s  cali¬ 
brated  using  10  years  of  monthly  data  on  ice  extent  .Sub¬ 
sequently,  the  satellite-derived  fields  of  ice  extent  h;iv(' 
been  used  by  Hibler  and  Walsh  (1982),  Hibler  and  .Ackley 
[1983],  Parkinson  and  Bindschadler  [1984],  Wai.sh  et  ai 
[1985],  Preiier  I19S51.  Hibler  and  Brj-an  11987).  Seminer 
[1987],  Stossel  et  al.  119891.  Fleming  and  Semtnor  (199!  1, 
and  Walsh  and  Chapman  (1991).  Satellite-derived  ice 
concentrations  in  the  vicinity  ofthe  Weddell  polyny  a  stimu¬ 
lated  the  modeling  studies  of  Martinson  et  a!.  (Ihsiy 
Parkinson  [1983],  and  Motoi  et  al.  !19871. 

While  most model-to-data comparisons oficc  extent  h.ive 
been  intended  to  assess  the  climatologie.s  simulated  by 
models,  several  studies  have  used  satellite-derived  data  to 
evaluate  the  ability  of  models  to  reproduce  mterannual 
fluctuations  of  ice  cover.  Walsh  et  al.  [1985],  for  exatnpic. 
found  that  the  Hibler  [1979!  two-level  model  reproduced 
approximately  50%  ofthe  year-to-year  vanance  ol  ire  cov¬ 
erage  in  20  longitudinal  sectors  of  the  .Arctic.  Wa!.-;h  and 
coworkers  used  intcrannually  varying  winds  and  air  tern 
peratures,  but  climatological  ocean  forcing  to  drive  a  .-^tand- 
alone  sea  ice  model.  Fleming  and  Semtnor  11991  j  subse¬ 
quently  showed  that  the  verification  slatistic.s  are  unproved 
somewhat  when  the  ice  model  i.s  coupled  to  an  interacDve 
ocean  model. 

More  specific  use  of  the  ice  concentrations  derived  froni 
passive  microwave  measurements  has  been  made  by  .Sin  izc 
et  al.  [19901,  who  used  an  ice  model  and  other  tool-  to 
diagnose  the  origins  of  a  large  area  of  reduced  icec  jncentr.)- 
tion.s  within  the  Arctic  pack  ice  dt  ring  1988.  .Serre.n  .ind 
coworkers  supplemented  the  passive  microwave  dat.i  w at  b. 
visible-band  imagery  during  the  sunlit  .vea.son. 
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Additional  uses  ot'satellite  measurements  for  model-to- 
data  comparisons  have  been  made  by  Walsh  and  Zwally 
i  19901,  whose  application  of  satellite-derived  multiyear  ice 
concentrations  is  discussed  in  Section  22.4.2,  and  by  Emery 
et  al.  [1991],  who  used  AVHRR  imagery  to  obtain  ice 
velocity  vectors  for  the  region  ea.st  of  Greenland.  In  the 
latter  study,  discrepancies  between  model-  and  satellite- 
derived  velocities  indicated  that  the  strength  of  an  ocean 
current  had  been  underspecified  in  the  model  forcing.  Ross 
and  Walsh  [1987]  compared  interannua!  fluctuations  cf 
Arctic  Ocean  surface  albedo  with  those  derived  from  satel¬ 
lite  imagery  by  Scharfen  et  al.  [19871. 

22.4.2  Hindcast  Studies:  Walsh  and  Zwally  [1990] 

Hindcaststudies  of  the  interannual  fluctuations  of  sea  ice 
have  been  performed  by  using  observational  analyses  of 
atmospheric  forcing  variables  iair  temperatures  and  pres¬ 
sure-derived  winds)  to  drive  the  Hibler  two-level  model  over 
time  scales  of  several  years  to  several  decades  [Walsh  et  al., 
1985;  Walsh  and  Zwally,  19901.  The  inclusion  of  multiyear 
ice  in  the  simulations  was  motivated  by  the  availability  of 
the  multiyear  ice  concentrations  derived  from  satellite 
passive  microwave  measurements.  The  oceanic  component 
of  the  mode!  consisted  of  only  a  constant-depth  mixed  layer 
and  time-independent  currents.  The  currents  were  ob¬ 
tained  geostrophically  from  a  prescribed  dynamic  topogra¬ 
phy. 

As  noted  in  the  previous  subsection,  runs  of  the  model  in 
this  mode  explain  about  50%  of  the  interannual  variability 
of  satellite-derived  coverage  of  total  sea  ice  (first-year  plus 
multiyear  ice).  The  simulated  values  of  multiyear  ice 
concentration  are  generally  larger  than  those  derived  from 
the  SMMR  data.  Although  several  arbitrary  elements  of 
the  model  formulation  may  contribute  to  this  disscrepancy, 
Comiso  [1990]  has  recently  reported  that  SMNIR-derived 
multiyear  concer.  trations  from  the  SMMRTeam  algorithm 
are  also  smaller  than  one  would  infer  from  the  concentra¬ 
tions  of  ice  surviving  the  summer  melt  in  the  Arctic,  The 
interannual  variations  of  SMMR-  and  model-derived 
multiyear  ice  coverage  generally  show  somewhat  ambigu¬ 
ous  agreement  in  the  Arctic  Ocean  (Figure  22-2). 

The  SMMR  data  also  permit  the  estimation  of  general¬ 
ized  resultant  velocity  vectors  on  the  basis  of  changes  in  the 
mu  Itiyear  ice  edge  over  periods  ofseveral  weeks  to  a  season. 
Such  estimations  are  limited  to  the  nonsummer  months, 
when  surface  wetness  does  not  contaminate  the  multiyear 
ice  signature. 

Instances  of  convergence/divergence  m  the  numerical 
model  results  have  been  confirmed  by  corresponding  changes 
in  the  multiyear  ice  concentrations  derived  from  the  SMMR 
data,  although  the  model-derived  changes  are  smaller  than 
those  obtained  from  SMMR.  Larger  samples  of  such  com¬ 
parisons,  including  buoy-derived  velocity  divergences,  are 
needed  to  determine  whether  there  is  a  bias  in  either  the 
SMMR  or  the  model  depictions  of  .•^hort-term  divergence 
events.  The  outcome  of  such  comparisons  has  important 


Fig.  2'2-2.  Timescne,sof  model-denvcd  multiyear  andrir^l-AeariC': 
coverage  in  the  Arctic.  Also  shown  are  N'ovembe.’'-April  averages 
from  SMMR  data,  adju.=ted  to  include  an  estimated  8  2x  1 C'  km-  m 
the  unobser\’ed  area  north  of  8-5°  N  [Walsh  and  Zwally,  lyooi 

implications  for  the  validity  of  the  simulated  mas.s  budget, 
since  most  of  the  ice  growth  in  the  central  i\rctic  occurs  in 
areas  of  thin  ice  or  open  water  (e.g.,  leads)  that  result  from 
divergence  of  pack  ice. 

22.4.3  The  Use  of  Models  With  Data  in  Forecasting 

The  Polar  Ice  Prediction  Syutem  is  a  sea-ice  forecasting 
system  based  on  the  Hibler  ice  model  [Preller,  1 985;  Prc-ller 
and  Posey,  1989].  It  is  used  for  the  daily  prediction  of  ice 
drift,  ice  thickness,  and  ice  concentration  in  the  Arctic  by 
the  U.S.  Navy  (Figure  22-3).  In  this  system,  the  Hibler  ice 
model  is  driven  b}'  daily  atmospheric  stresses  and  heat 
fluxes  from  the  Navy  Operational  Global  Atmospheric  Pre¬ 
diction  System  and  by  monthly  mean  ocean  currents  and 
heat  fluxes  derived  from  the  Hibler  and  Bryan  [  1 98T  j  icc-- 
ocean  model.  The  model  is  run  in  an  operational  environ¬ 
ment  at  the  Fleet  Numerical  Oceanography  Center '  P'N OC  ' 
A  120-hour  forecast  of  several  different  ice  fields  is  produced 
daily.  Each  day’s  forecast  is  initialized  from  the  pre\  iou- 
day’s  24-hour  forecast  (Figure  22-4).  Once  each  week,  the 
model’s  ice  concentration  field  is  initialized  from  a  digitized 
field  of  observed  ice  concentration  made  available  by  the 
Navy/NOAA  Joint  Ice  Center  (Figure  22-5).  This  field  ,t 
subjective  analysis  derived  from  multiple  data  .sourri  s 
fWohl,  1991  ].  Satellite  data  u.sed  in  the  analysis  come  from 
the  AVHRR,  visible  data  from  the  DMSP  OL.S.  and  pas.sivc 
microwave  data  from  the  SSM/I.  In  addition  to  these 
satellite  data,  any  available  ice  reconnaissance  fiighl.'  .a 
ship  observations  arc  blended  into  the  weekly  analA  si.-. 

During  model  initialization,  the  ob.serc-ed  ice  coiu  cntra 
tion  totally  rcplace.s  the  modeled  ice  concentration  from  tin 
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previous  day’s  24-hour  forecast.  Other  variables  from  the 
previous  day’s  forecast  remain  unchanged,  with  the  excep¬ 
tion  of  the  ice  thickness  and  the  heat  stored  in  the  ocean’s 
mixed  layer.  These  fields  are  adjusted  only  near  the  ice  edge 
to  be  consistent  with  the  ice  edge  defined  by  the  JIC  data  as 
shown  in  Figures  22-6(a)  and(b).  This  adjustment  involves 
removing  ice  from  regions  that  should  not  be  ice  covered  and 
adding  heat  to  the  mixed  layer  by  raising  its  temperature. 
If  the  model  forecasts  no  ice  in  a  region  where  the  data 
indicate  it  should  exist,  heat  is  removed  from  the  mixed 
layer  by  setting  the  temperature  back  to  freezing  and  either 
0.5  or  1 .0  m  of  ice  is  added  to  that  grid  cell  depending  on  the 
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Fig.  22-6.  PIPS  24-hour  foreca.st  ice  concentration  from  (a)  June  G, 
1991,  and  (b)  June  7, 1 991 ,  immediately  following  the  initialization 
of  PIPS  by  the  JIC  analysis.  (The  USSR  is  now  referred  to  as  the 
Commonwealth  of  Independent  States.) 


observed  concentration  [Preller  and  Fo.sey,  1989]  Two 
similarly  designed  forecast  systems,  the  Regional  Polar  Ice 
Prediction  System  Barents Sea<RPIPS-B]  and  the  Regional 
Polar  Ice  Prediction  System  Greenland  Sea  (RIMPS-G* 
forecast  over  the  Barents  and  Greenland  Sea;,  at  a  higher 
resolution  than  PIPS  [Preller  ct  al.,  1989] 

Data  assimilation  and  model  initialization  in  an  opera¬ 
tional  forecasting  system  are  dependent  on  a  number  of 
different  factors.  Availability  of  the  data  is  the  first  impor¬ 
tant  factor.  The  data  must  be  available  Ui  the  computer 
system  used  by  the  forecast  model. 

The  second  factor  is  the  age  of  the  data.  The  model  uses 
data  to  provide  the  most  realistic  initial  state  for  the 
forecast.  If  the  data  are  old,  they  can  provide  incorrect 
information  to  the  model.  The  JIC  completes  its  anaiysns 
each  Thursday;  the  analysis  is  ready  for  use  in  PIPS  by  th-,.- 
end  of  the  day  on  Friday.  The  two-day  lag  in  JIC  data  is  not 
optimal,  but  is  considered  acceptable  in  the  PIPS  model. 
PIPS  uses  a  resolution  of  127  km.  The  ice  edge,  on  average, 
does  not  move  more  than  100  km  over  48  hours,  so  the 
observ  ed  changes  may  not  be  resolved  by  the  model.  There 
are  certainly  exceptions,  such  as  the  appearance  of  the 
Odden  over  a  24-hour  period  or  the  rapid  movement  of  the 
ice  edge  brought  about  by  extremes  in  the  atmospheric 
forcing(strongwinds).  A48  hourlimedelayin  thedafacan. 
however,  be  a  problem  in  the  RPIPS-B  and  RPiPS-G  mod¬ 
els.  Since  these  models  use  resolutions  of  20  and  25  km. 
respectively,  movement  of  the  ice  edge  in  these  regions  can 
be  a  grid  distance  or  moreover  a  48-hour  period  Ideally,  the 
data  used  for  initialization  should  be  no  more  than  24  hours 
old.  The  importance  of  the  age  of  the  data  alsoimplie.s  that 
data  assimilation  techniques  such  as  Kalman  filiering. 
which  fits  a  time  series  of  data  into  a  predicti%e  model,  may 
not  provide  the  most  recent  ice  conditions  for  initializing  a 
forecast.  Also,  the  Kalman  filter  may  not  be  the  optimal 
choice  in  an  operational  system  where  computer  time  i.s 
limited  (  the  Kalman  filter  is  computer  intensive )  and  where 
historical  data  arc  not  often  saved  and  available  for  u.se  with 
the  filter. 

The  third  important  factor  is  the  frequency  of  the  data. 
The  JIC  analysis  is  available  only  once  per  week.  How'over, 
FNOC  now  has  a  digitized  data  set  of  ice  concent  ration  from 
the  SSM/1  available  in  real  time  each  day  (Figure  22-7). 
That  data  set  is  presently  being  tested  as  an  initialization 
field  for  PIPS  on  those  Jay.s  that  the  JIC  analysis  is 
unavailable. 

In  addition,  Cheng  and  Preller  [por,sonal  communication, 
1991)  have  designed  a  method  of  blending  two  or  more 
concentration  fields  into  one  maximum  likelihood  estimate 
using  nonlinear  regression.  The  technique  can  be  described 
by  borrowing  the  terminology  of  the  physical  and  measure¬ 
ment  models  from  Thorndike  119881  anj  Thomas  and 
Rothrock  11989)  (see  Chapter  23).  The  physical  mode! 
de.scribing  pack  ice  evolution  can  be  written  as 

A,  +  o  (9. 
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Fig.  22-7.  SSM/I  ice  concentration  from  FNOC,  derived  using  the 
Navy  algorithm,  for  June  7, 1991.  Values  have  been  interpolated 
to  and  plotted  on  the  PIPS  grid. 


where  X  is  the  ice  concentration,  H  is  the  physical  model 
PIPS,  e  is  the  residual  of  the  physical  model,  i  and  j  are  the 
spatial  indices,  and  t  is  the  time  index.  The  physical  model 
uses  dynamic  and  thermodynamic  equations_^long  with 
atmospheric  and  oceanic  forcing  to  predict  Xj‘*^.  The 
measurement  model  is  used  to  update  the  value  ofZ  for  the 
physical  model.  Written  in  a  general  form, 

Zfn  =  Gm{Xn^  (10) 

where  Z  represents  the  data  (i.e.,  observations,  measure¬ 
ments,  constraints,  and  assumptions),  G  is  the  measure¬ 
ment  model  (in  this  case  the  continuity  equation  for  ice 
concentration),  e  is  the  residual  between  the  measurement 
and  the  corresponding  model  value,  and  m  and  n  are  spatial 
indices. 

The  JIC  analysis  and  passive  microwavedata(otherdata 
may  be  included)  provide  field  observations  or  measure¬ 
ments.  The  PIPS  ice  concentration  provides  a  data  set  of 
historical  information,  while  the  continuity  equation  for  ice 
concentration  represents  the  constraint  for  these  ice  con¬ 
centration  fields.  All  of  these  data  and  the  constraints 
imply  that  the  measurement  model  is  defined  as  an  over¬ 
determined  problem,  which  is  particularly  suited  to 
nonlinear  regression.  The  method  also  allows  for  emphasiz¬ 
ing  the  reliability  of  the  data  in  regression  by  applying 
variable  weights  to  the  data. 


Remotely  sensed  data  have  also  been  used  for  validating 
the  Navy’s  operational  sea-ice  forecast  systems.  Each  IJ.S. 
Navy  operational  model  is  put  through  an  operational 
evaluation,  lasting  anywhere  from  three  months  to  a  year. 
Since  1987,  an  operational  forecast  of  at  least  24  hours  has 
been  made  daily  and  archived  by  the  Naval  Research 
Laboratory  (NRL)  (formerly  the  Naval  Oceanographic  and 
Atmospheric  Research  Laboratory).  These  forecasts  are 
identical  to  those  made  by  h’NOC  and  are  used  by  NRL  for 
continuing  model  verification.  During  the  operational 
evaluation  of  PIPS,  the  model-derived  ice  edge  was  com¬ 
pared  qualitatively  and  quantitatively  to  the  ice  edge  de¬ 
rived  by  JIC.  This  study  showed  that  when  initialized  at 
least  once  each  week,  the  mean  error  in  the  forecast  i'e  edge 
location,  after  seven  days,  was  less  than  or  equal  to  one  grid 
distance.  This  error  increased  if  the  initialization  took 
place  less  frequently  IPreller  and  Posey,  1989]. 

In  addition  to  ice  edge  verification,  the  PIPS-derived  ice 
drift  is  verified  against  Arctic  buoy  data  on  a  regular  basis. 
Ice  drift  dertved  from  the  change  in  buoy  location  is  statis¬ 
tically  compared  to  PIPS  ice  drift.  The  result  of  this 
verification  has  been  the  continued  improvement  in  the 
accuracy  of  the  PIPS  ice  drift  field,  mainly  due  to  improve¬ 
ments  in  the  treatment  of  surface  wind  stress  on  the  ice 
[Preller  and  Posey,  1989].  During  the  operational  evalua¬ 
tion  ofboth  regional  models,  satellite-derived  ice  drift  data 
from  successive  AVHRR  data  were  used  for  the  validation 
ofthe  model-derived  ice  drift.  Figures 22-8<a)  and  (b).  In  the 
case  ofboth  regional  models,  the  direction  ofthe  modeled  ice 
drift  agreed  well  with  the  satellite-derived  ice  drift,  how¬ 
ever,  the  magnitude  of  the  modeled  drift  was  weaker  than 
that  derived  from  the  satellite  data. 

22.5  StwiARv  AND  Conclusions 

Numerical  models  and  observational  data  have  been 
used  together  in  a  number  of  ways.  Data  obtained  from 
drifting  ice  station,  aircraft,  and  ship  measurements  col¬ 
lected  over  many  years  provided  highly  localized  point 
measurements  that  were  used  to  define  many  characteris¬ 
tics  of  sea  ice.  A  more  synoptic  view  of  sea  ice  was  obtained 
from  the  data  collected  in  field  experiments  such  as  AIDJEX 
and  MIZEX.  These  data  provided  the  basis  for  an  improved 
understanding  of  the  dynamic  and  thermodynamic  pro¬ 
cesses  that  affect  sea  ice.  As  a  result,  more  sophisticated 
formulations  for  ice  motion,  including  the  effect  cf  internal 
ice  stress  and  the  growth  and  decay  of  sea  ice,  have  been 
determined  and  incorporated  into  numerical  models.  The 
models,  in  turn,  are  used  to  help  provide  insight  into  the 
accuracy  of  these  formulations  by  validating  the  model 
results  against  observations. 

Satellites  have  been  able  to  provide  a  source  of  observa¬ 
tional  data  both  over  larger  areas  and  for  longer  periods  of 
time  than  the  more  conventional  field  ob.servations  A 
combination  of  fine-resolution  satellite  imagery,  such  as 
visible  or  infrared  imagery,  and  radar  altimetry,  along  with 


coarser  resolution  data  such  as  that  from  the  passive 
microwave  sources,  have  been  used  for  verifying  and  vali¬ 
dating  many  numerical  models.  A  comparison  of  the  ice 
edge  derived  frozn  these  sources  vith  the  modeled  ice  edge 
has  been  the  favored  method  of  model  validation. 

Satellite  passive  microwave  data  are  particularly  appeal¬ 
ing  to  the  large-scale  models  used  to  study  interannual 

(a) 


Fig.  22-8.  Ico  drift  derived  from  (al  successive  AVFIKRimages  from 
March  2'i  and  24.  1  989,  and  (l)t  the  24 -hour  forecast  of  RPIPS-B  for 
the  same  time  period. 


variability  of  .sea  ice.  Passive-riucrowave  total  ice  coticeii- 
tration  has  been  used  in  the  majority  of  model-io-dai.i 
comparison  studies.  However,  avaihibility  of  niultiyear  ice 
concentration  from  the  passive  microwave  data  motivated 
Walsh  and  Zwalty  II 9901  to  incorporate  multiyear  lee  into 
a  version  of  the  Hibler  ice  model.  Using  thi.s  formulation, 
they  were  able  to  .show  some  conipaiibilily  between  the 
interannual  variations  of  the  model-  and  SMMK-derived 
multiyear  ice  concentration. 

A  number  of  sources  of  remotely  .sensed  data  are  m-sed  to 
create  a  weekly,  subjective  analysis  of  ice  concentration  by 
JIC.  This  analysis  is  used  as  a  weekly  initialization  field  for 
each  of  the  three  sea-ice  forecast  .systems  run  by  the  U  S 
Navy;  PIPS,  RPIPS-B,  and  RPIPS-G,  Two  upgrades  are 
scheduled  for  the  Navy  forecast  systems  in  the  near  future. 
The  first  upgrade  involves  the  use  of  passive  microwave 
data  as  an  initialization  field  for  the  forecast  systom.s 
Passive  microwave  data,  now  available  in  real  time  from 
FNOC,  are  being  tested  as  a  data  source  that  could  provide 
more  frequent  initialization  of  these  forecast  systems.  The 
second  upgrade  involves  the  replacement  of  PIPS  with  the 
PIPS  2.0  model.  PIPS  2.0  will  provide  ice  forecasts  for  mo.'i 
of  the  ice-covered  regions  in  the  Northern  Hemisphere  on  a 
half-degree  grid  (Figure  22-9).  The  design  of  this  .system 
will  be  similar  to  that  of  PIPS,  but  will  include  a  fully 
coupled  ice-ocean  model.  Along  with  the  JIC  analysis,  the 
SSM/I  ice  concentration  will  provide  real-time  initializa¬ 
tion  data  over  the  entire  model  domain. 

In  addition  to  the  passive  microwave  data,  the  recent 
launch  of  the  ERS-1  satellite  is  now  providing  finc-resolu- 
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F'ig.  22-9.  FIP.S  2.0  moiithlv  mean  Ice  concentration  for  M,iri;!i 
1986, 
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tion  coverage  of  sea  ice  conditions  via  the  Synthetic  Aper¬ 
ture  Radar.  Considerable  potential  exists  for  the  use  of  SAR 
data  in  model  studies  of  the  parameterization  of  thin  ice  and 
open  water  regions.  In  addition,  SAR  may  help  provide 
more  realistic  parameterizations  of  the  surface-atmosphere 
heat  exchange,  thereby  serving  as  a  stimulus  to  the  cou¬ 
pling  of  ice  and  atmosphere  models,  as  well  as  atmosphere— 
ice-ocean  models. 

Although  large  quantities  of  remotely  sensed  data  al¬ 
ready  exist,  there  is  great  expectation  for  what  will  be 
available  in  the  future.  It  is  anticipated  that  remotely 
sensed  data  will  provide  additional  important  information 
to  models  on  such  variables  as  albedo,  ice  surface  tempera¬ 
ture,  snow'  cover,  ocean  surface  temperature,  and  atmo¬ 
spheric  properties  such  as  cloud  cover,  wind  speed,  and 
radiation  budgets. 

Although  the  interaction  between  data  and  models  has 
been  ongoing  for  many  years,  the  future  would  suggest  an 
even  greater  opportunity  for  interplay  between  the  two. 
The  combination  of  new  sources  of  remotely  sensed  data 
and  new  products  derived  from  these  data  should  assist  in 
the  evolution  of  improved  formulations  for  dynamic  and 
thermodynamic  processes  affecting  sea  ice.  These  new 
formulations  will  be  used  to  create  more  sophisticated 
models  of  the  interaction  between  atmosphere,  ice,  and 
ocean.  In  addition,  recent  work  by  Thomas  and  Rothrock 
[1989],  in  Chapter  23,  and  the  work  presently  being  done  by 
the  U.S,  Navy,  suggests  that  the  assimilation  of  data  into 
models  can  provide  both  improved  models  as  well  as  im¬ 
proved  forecasts  of  sea  ice  conditions. 
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